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Экстраполяционные блочные 
одношаговые численные методы 
решения жестких задач Коши 


Экстраполяционные методы решения задачи Коши ДлЯ обыкновенных дифференциальных уравнений 
обладают высокой степенью потенциального параллелизма. Данная статья посвящена разработке и 
анализу эффективности параллельных алгоритмов локальной экстраполяции на базе явных опорных 
методов. Разработанные алгоритмы реализованы на параллельных системах с распределенной памятью 
и топологией гиперкуб. Получены оценки времени выполнения и обменов, общих накладных расходов 
на параллелизм, ускорения и эффективности параллельного решения. 


Введение 


Экстраполяционные методы Ричардсона предназначены для получения высоко- 
точного решения задачи Коши и интегрирования систем обыкновенных дифферен- 
циальных уравнений (СОДУ) со сложными правыми частями. Практическое приме- 
нение экстраполяционных методов затруднено в связи с большой вычислительной 
сложностью их последовательных реализаций. Кроме того, использование техноло- 
гии локальной экстраполяции на основе явных опорных схем ограничивает исполь- 
зование методов областью нежестких задач. Поэтому построение эффективных 
параллельных блочных неявных методов локальной экстраполяции Ричардсона — 
один из наиболее реальных способов сокращения времени интегрирования много- 
мерных жестких начальных задач. 

Целью данной работы является разработка параллельных блочных экстра- 
поляционных одношаговых методов численного решения жестких задач Коши и исследо- 
вание эффективности их отображения на параллельные вычислительные структуры с 
распределенной памятью. 


Общая характеристика 
технологии локальной экстраполяции 


Решение задачи Коши для систем обыкновенных дифференциальных уравне- 
ний первого порядка: 


ах вр. 
ты = Л(х,У(х)), (1) 
У(хо)= Уд, 


рассматривается при переходе из точки х, в точку х„/=х,+Н, где Н- базовая 
длина шага (рис. 1). Выбирается ряд натуральных чисел Р = {п,п.,..., Пу, ...}, 


п] < П> <...< Пу <... и, соответственно, последовательность шагов интегрирования: 
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п >12 >...> М> >... Где р; =Н/п;. Задается опорный численный метод 
порядка и) и вычисляются приближенные решения исходной задачи Коши в точке 


Х1+1: 1 = Ув (Хи +Н)=УИИ=ИК. 


х х +Н 


п п 


Рисунок 1 — Схема технологии локальной экстраполяции Ричардсона 


Выполнив вычисления для ряда последовательных значений 1, по рекуррент- 
ному соотношению определяют значения для произвольных 1/ по схеме локальной 


полиномиальной экстраполяции Эйткена-Невилла: 


Т. . Е а о 2Т-ы 
+1 1] (п, п,” =] 


Величина БР равна 1 в общем случае, для симметричных опорных методов В 


равно двум. Т, — приближенное решение задачи Коши, полученное численным 


методом порядка -+(/-1)-Б с шагом й;. Достоинство этого метода состоит в том, 
что он дает целую таблицу результатов вычислений (табл. 1), которые образуют 
последовательность вложенных методов, позволяют оценить локальную погреш- 
ность и выбрать стратегию для управления шагом интегрирования [1-3]. 


Таблица 1 — Экстраполяционная таблица 


1) +В ... п) +(-ПБ 
Ти 

То Т›2 

Те Тр ... Ти 


Построим эффективный с точки зрения минимизации вычислительных затрат 
последовательный метод с использованием технологии Ричардсона. Так как экстра- 
поляционная технология включает: опорный численный метод решения задачи Коши, 
последовательность сеток интегрирования, рекуррентное правило вычисления значе- 
ний приближенного решения, то эффективность ее применения напрямую зависит от 
правильного выбора и сочетания всех трех составляющих. 
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Для получения множества шагов интегрирования применяются числовые по- 
следовательности, образованные гармоническим рядом, рядом четных чисел, степенями 


двойки [4-6]: 1) Р, = { 1,2,3,4,5,6,..}; 2) Р, ={1,2,4,6,8,10,12,...1; 3) Р, = {1,2,4,8, 
16,32, 64,...}.Вычислительная сложность использования последовательностей нату- 
ральных чисел для произвольных опорных методов определяется числом обращений 
к функции / - правой части системы исходных ОДУ: 


М(К) = п, + п, +...+п,, 
РР: МЕ) =(К? +К)/2, 
Р: МЕ) =Е?-К+1, 
ВР: М®=2-1, 


где А — число строк экстраполяционной таблицы. 
При одном и том же опорном методе наиболее затратной является последо- 
вательность Р;, которая фактически повторяет процесс удвоения шага; последова- 


тельности В; И Р, имеют приблизительно одинаковую вычислительную сложность. 


Выбор опорного метода интегрирования базируется на исследовании влияния 
порядка метода и его свойств на объем необходимых вычислений [4-7]. Следо- 
вательно, если вычислены А строк экстраполяционной таблицы, то для симметрич- 
ных опорных методов имеем Г,, в качестве аппроксимации наивысшего порядка 


точности, равного 2%, и, кроме того, точность аппроксимации Т,_/ ;‚ имеет порядок — 
(2К—2). Соответственно, для произвольных опорных методов имеем: Ти — (К-Г)и 


Т.к - (К -2) порядка точности. Для определения апостериорной локальной погреш- 
ности и управления шагом интегрирования используется величина: |7 ТК ТЕК | . Сум- 


мируя результаты, можно сделать вывод, что для уменьшения накладных вычисли- 
тельных расходов при применении технологии локальной экстраполяции следует 
выбирать опорный метод малого порядка с достаточными свойствами устойчивости. 

Особенностью интегрирования задачи Коши для СОДУ по технологии ЛЭР 
является тот факт, что экстраполяционная таблица вычисляется для каждой размер- 
ности исходной системы. Теоретически нет никаких ограничений на длину экстрапо- 
ляционной таблицы, но конечность памяти машины и накапливание ошибок округления 
ограничивают длину экстраполяционной таблицы сверху и, как правило, используется: 


0-20 


К тах < 10. Поэтому для получения высокоточных приложений (1 Иа ), где вы- 


тах — 

числения правой части СОДУ достаточно сложны, не имеется другой возможности, 
как использовать технологию локальной экстраполяции Ричардсона на базе опорного 
метода высокого порядка точности. 

Потенциально вычисления по технологии локальной экстраполяции содержат 
три источника внутреннего параллелизма: 

— системный параллелизм (ограничен размерностью СОДУ, т); 

— параллелизм экстраполяции (ограничен размерностью таблицы экстрапо- 
ляции, К); 

— параллелизм опорного метода (малая степень параллелизма). 
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Технология локальной экстраполяции на основе 
одношаговых блочных неявных методов 


Блочные, или многоточечные, методы решения задачи Коши особенно акту- 
альны, поскольку хорошо согласуются с архитектурой параллельных вычислитель- 
ных систем и не требуют вычисления значений в промежуточных точках, что значи- 
тельно повышает эффективность счета. В то же время одношаговые блочные методы 
обладают высокой устойчивостью, что позволяет применять их для решения жест- 
ких задач Коши [8], [9], так как позволяют получать решение одновременно в 
нескольких точках сетки интегрирования. 


Множество точек равномерной сетки 42, : { х; } ]=1,М разбивается на № 
блоков. Каждый блок содержит К точек и при этом М < М. Предполагается, что в пре- 
делах блока все точки равноудалены друг от друга: 


ХЕХ + =ЬК, (2) 


где /— номер точки в блоке, 1=/,К; п- номер блока, п=/,М№; х„,- точка с номе- 
ром Г, принадлежащая блоку и; х„)-— начальная точка и-го блока; Х„; - конечная 
точка и-го блока. 
К 
Множество точек и-го блока из К точек обозначается Т““”. При этом имеет 


место равенство: Х„; =Х,„. 10. Пусть у.о — приближенное значение решения задачи 


Коши в точке Хх, ) — начальной точке обрабатываемого блока. 


Разностные уравнения для блока из А точек могут быть записаны в виде: 
[9 ее ЗЫ 
Упя = Уп,0 + Ед + У ан; = РК; п = М . (3) 
7 

Разложением в ряд Тейлора входящих в невязку функций можно показать, что 
одношаговый „-точечный блочный метод имеет наивысший порядок аппроксимации, 
равный АК+/, следовательно, локальная ошибка в узлах блока имеет порядок 
О( ре ). Блочные параллельные методы относятся к классу неявных, поэтому для 


вычисления приближенных значений решения задачи Коши необходимо разрешить 
систему нелинейных уравнений. Одним из способов получения решения является 
метод простой функциональной итерации: 


УЕ +ЙЕ, ЕЁ , п=12,..,М, 


. —_ (4) 
Ут > Ут,0 о И(Б,Е о + рр. > 1 == 0, Г а р, 


= 
где и — номер блока, и = /,2,...,№М; 1— номер точки блока, #1= /1,К; [ — номер текущей 
итерации /[=0,2-1; [Г - максимальное число ненулевых итераций. 


Реализация технологии локальной экстраполяции Ричардсона для блочных 
методов требует многократных вычислений на одном и том же интервале интег- 


рирования с использованием опорного блочного Ку-точечного метода порядка 7, на 
сгущающихся равномерных сетках: 


а) Зи, [хр л=ЬМ, с шагом й; =Й/п; в М, блоках; 
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6) 2, Ех =тм, с шагом й› =Й/п> в №, блоках; 

в) (т, У с шагом Й;, =Й/и; в №, блоках, где К-—число строк 
экстраполяционной таблицы. 

Базовый шаг интегрирования равен Й, =й, и, = [, то есть основой счета является 
сетка: {2,. Для генерации вспомогательных сеток выбирается гармонический ряд 
Р/={пП›,...,Пу,..4=$2,3,4,5,...}, как наименее затратный в случае опорного метода 


произвольного типа. Блочный опорный метод должен иметь малый порядок точности: 


Ко — 
Тир = Ув; = Упо + БЕ ,0 + Хау 1 = То, Ко $4, (5) 
11 
так как с ростом И) вычислительные затраты на технологию в целом существенно 


возрастают, несмотря на линейное уменьшение длины экстраполяционной таблицы. 
Для неявных методов, к которым относятся и рассматриваемые блочные методы, это 
особенно важно, поскольку увеличение порядка, а следовательно, и числа точек бло- 
ка, влечет за собой увеличение порядка многократно решаемых СНАУ. 

Схема разбиения базового интервала на блоки при применении технологии 
локальной экстраполяции для двухточечного опорного метода и гармонической пос- 
ледовательности приведена на рис. 2. Заметим, что для блочных методов базовый интер- 
вал интегрирования — это некоторый блок с номером и основной сетки длины Н = Кой. 


Общее число точек, вычисляемое с разными шагами интегрирования й; =й/1,1 = 1,К и 
определяющее узлы сеток 42, по Р,, равно М; =Н/Й, = (Кой) /Й; = КоЁ. Соот- 
ветственно, количество блоков 1-й сетки №; = М; / № = ип; =М,1=1,К. 


Поскольку экстраполируются значения в узлах базовой сетки, необходимо ус- 
тановить механизм соответствия узлов сетки (2, основного счета узлам сеток для 


экстраполяции ке =2.. Для двухточечного опорного блочного метода и гар- 
монического ряда имеем: 

[) 1- четное число: х,.(й) <=> Х12 (1); 

2) 1- нечетное число: х,.(й) <> хай 

3) для любых 1: х„2(й) => х,12(й;). 


Значения первого столбца экстраполяционной таблицы определяются на 
основе следующих формул: 


К И 
ое Е , 
ТИ ы У’ ы У + - БЕ 0 т ани = й 0; П — р, ув 
, ) Е 
(6) 
®Ф_ ® и и а аа 
Ти У г РЕ: РЕ, т >.а, Ел 1 = [, 0; И = /, К 
= 
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и х 
р п,1 р п,2 
Хи,0 ——»—»—»<_———— МЫ -— 52, 
х % № х 
2 п,1 п,2 п+1,1 п+1,2 
Х,0 е—_—_—_— 6—8 - 22 
й/3 
Хи, 1 Хи, Хт+1,1 Хи+1,2 Хи+2,1 Хи+2,2 
Х,0 ее ——_— 2 —_— 6—6 —— 6 КФ 
Хи Хн,2 Хи+ 11 И Хин Ан 


Рисунок 2 — Схема разбиения на блоки по технологии локальной экстраполяции для 
одношаговых блочных методов 


Каждая ИЗ аппроксимаций решения получается за счет № ; раз примененной 


схемы одношагового блочного К -точечного метода с разными шагами интегриро- 


вания: 


т, в Г = М, у р (п), 


5,1 (в) 
(7) 
= М, -Т” (ВК), 


70 о. 
где Г, (й/1),1 =1,К — время, необходимое для решения задачи Коши для ОДУ опор- 
ным методом порядка И, с шагом Й,. Затем по формуле Эйткена-Невилла вычис- 
ляются приближения Т,, ‚Г..,... и Тур. 


Построение параллельных алгоритмов решения нелинейной задачи Коши для 
СОДУ по технологии Ричардсона базируется на иерархической декомпозиционной 
методике. После того как определены все составляющие технологии экстраполяции, 
процесс решения можно разбить на две последовательно выполняемые подзадачи: 

1) вычисление А аппроксимаций решения в точке х„./ с разными шагами ин- 


а 


тегрирования: У», (х„+Н)=У, 4 |; 

2) построение экстраполяционной таблицы. 

Организация параллельных вычислений может производиться двумя спосо- 
бами, отсюда и две различные макрооперационные схемы алгоритма. Первый вари- 
ант подразумевает использование только системного параллелизма, второй — комби- 
нацию параллелизма экстраполяции и системного. Рассмотрим первую подзадачу 
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алгоритма. В качестве основной макрооперации для обоих вариантов вводится — 
вычисление аппроксимации точного решения в точке сетки с заданным шагом интег- 
рирования на базе явного опорного метода. Граф влияния для первого варианта 
макрооперационной вычислительной схемы приведен на рис. 3. 


(1) (Е-1) (К) 


Уз — АН рт. Уна = 


(Е-1) (К) 


(1) 
Ут ИВ Ут — Ут —_ 


й1 пу 


Рисунок 3 — Граф влияния первой макрооперационной схемы ЛЭР 


(системный параллелизм) 


При таком варианте все аппроксимации решения вычисляются последователь- 


но одновременно для каждой размерности системы. По второй макрооперационной 
схеме параллельно вычисляется К независимых аппроксимаций решения в точке 


х, + Н , граф влияния приведен на рис. 4. Второй вариант алгоритма имеет большую 


степень параллелизма, равную Рор=т-К, в отличие от первого, для которого: 


Рор =т. Однако этот недостаток схемы № 1, а именно последовательное выполне- 


ние макроопераций, может быть компенсирован лучшей балансировкой загрузки 
многопроцессорной вычислительной системы, а также низкой интенсивностью обменов. 


122 


55(1) 


У 


(2) 


о о,’ 


=7(3) 
Ул+1 


й, 5) 
Ф_>е_>е_> ”,-! 


Рисунок 4 — Граф влияния второй макрооперационной схемы ЛЭР 
(параллелизм экстраполяции) 
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Вторая подзадача, а именно вычисление экстраполяционной таблицы, в зависи- 
мости от способа решения первой подзадачи, также может быть распараллелена 
разными способами. В качестве макрооперации для этой подзадачи вводится вычис- 


ление экстраполированного значения 7 „/=1т1=2,К;|=Ьк, вычисленного по фор- 


муле Эйткена-Невилла (рис. 5). 


РР ай — 
7 т =’ Ве 
Г т. = 
Ту ть РЯ 
Ти ее”. 


Рисунок 5 — Вычисления по схеме полиномиальной экстраполяции 
Эйткена-Невилла, и =4 


Распараллеливание построения таблицы экстраполяции может быть осуществ- 
лено за счет независимого выполнения вычислений по формуле полиномиальной 
экстраполяции для каждой размерности системы (Дор=т ), а также за счет рас- 


пределенного вычисления каждой строки таблицы (Рор=^А-1). Как правило, 


т>к › Поэтому появляется возможность совместить оба вида параллелизма. На рас- 
смотренном этапе детализации дан предварительный анализ эффективности каждой 
из вычислительных схем, перейдем к использованию внутреннего параллелизма метода. 

Первая макрооперационная схема алгоритма локальной экстраполяции практи- 
чески повторяет макрооперационную схему для правила Рунге, с тем различием, что 
теперь у нас не три, а К последовательно вычисляемых аппроксимаций по ЯМРК. 
Заметим, что каждую макрооперацию схемы локальной экстраполяции можно пред- 
ставить как п; раз используемую макрооперацию однократного применения ЯВНОЙ 


численной схемы одношагового метода. 

Далее по декомпозиционной технологии необходимо распределить данные по 
процессорам. Если размерность процессорного поля велика и больше, чем длина экс- 
траполяционной таблицы, появляется возможность выполнять макрооперации первой 
подзадачи параллельно. Для этого необходимо разбить процессоры на группы таким 
образом, чтобы обеспечить наиболее рациональную балансировку загрузки много- 
процессорной системы. Рассмотрим следующие способы разбиения процессоров на 
группы, а соответственно и распределения данных по процессорам [4-7]: 1) равно- 
мерный, 2) пропорциональный, 3) комбинационный. 

При наиболее простом, равномерном способе каждый процессор или процессор- 
ная группа будут отвечать за вычисление равного количества аппроксимаций. Если 
размерность процессорного поля равна р, количество строк в экстраполяционной 


таблице — А, то при таком способе будет А групп по р, = Гр/ К | 1 = 1,К процес- 
соров в каждой и каждый процессор в 1-й группе будет содержать т, =|т/ р, | 
компонент соответствующей аппроксимации решения и соответствующего экстрапо- 
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лированного значения. Для равномерного разбиения: р, =р, =...=р, и, соответс- 


твенно, Ш, =т, =... =т,. Каждая группа будет отвечать за вычисление определен- 


К 

ной аппроксимации решения: Г,,, / = .й, время решения определяется по тах Т;у, 
= 

и все группы, кроме последней, будут простаивать некоторое время. 

Рассмотрим различные варианты организации параллельных вычислений экс- 
траполяционной таблицы по схеме Эйткена-Невилла при равномерном распределе- 
нии для вычислительной схемы № 2. После параллельного вычисления К аппрокси- 
маций решения на А группах процессоров, имеется следующее распределение исход- 
ных данных для вычисления экстраполированных значений. Первая группа процес- 


соров содержит вектор Т,,,;/ = 1, т, вторая — Т,,,;/=1,т и, соответственно, К -я — 


Т р] И 7 =/ ‚т, каждый процессор в группах содержит одинаковое т; число компонент. 


Для вычислений элементов экстраполяционной таблицы могут быть реали- 
зованы два варианта. Первый состоит в использовании системного параллелизма, 
каждый процессор группы передает соответствующему процессору соседней группы 
все компоненты своего вектора аппроксимации решения. Так, первый процессор первой 


группы передаст первому процессору второй группы подвектор: Т,,,,/ = ‚т, для 


вычисления такого же подвектора экстраполированного значения: Т,, ‚, / = 1, т, ‚ второй 
процессор первой группы передает Т,,,, / = т, + 1,2т, второму процессору второй 
группы для вычисления Т,,,,/ = т, + [,2т, и так далее, первый процессор (К -—1) 


группы передает первому процессору А группы значения подвектора Т,/,,/= 


=т,. для вычисления Т,,,,/=1,т,. Высота такого параллельного алгоритма 


равна: А-—/. То есть на первом шаге вычисляются подвектора второй строки таб- 
лицы экстраполяции, на последнем А-—1 шаге вычисляются соответствующие под- 
вектора К -й строки таблицы экстраполяции. 

Второй алгоритм вычисления экстраполяционной таблицы состоит в следую- 
щем: в каждой группе процессоров ровно один процессор будет отвечать за передачу 
данных между группами. Для этого после подсчета элементов первого столбца 
таблицы в каждой группе процессоров производится обмен значениями по типу 
«все-всем», в результате каждый процессор в группе будет содержать весь вектор 


аппроксимации решения, а не его часть. Процессоры первой группы -— Т,,, второй — Г), 
и ^-й -— Т,,. Затем 1-й процессор каждой группы, кроме последней, передает всем 


процессорам соседней группы вектор решений для подсчета экстраполированного 
значения. Далее каждый процессор каждой группы, кроме первой (ширина алгорит- 
ма уменьшается с каждым шагом), считает свой подвектор (длиной т, } своего вектора 


экстраполированного значения: процессоры второй группы -— Х,, , третьей — Т., и К-Й-— 
Т,›. На последнем А -/ шаге будут работать только процессоры последней группы, и 


каждый будет считать т, значений вектора Ти, . Чтобы увеличить балансировку загруз- 
ки, разобьем процессорное поле на такое же число групп, но не равномерно, а пропор- 


м [Я 
ционально числам и,,1= 1: ру: р>:...:рьЕ=т:т:.гт, Ур;=р, р; т .п,. Алго- 
‚= 
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ритм разбиения на группы следующий, если число процессоров не кратно №(А), то 


Р 
МК) 
каждую группу, начиная с последней. Таким образом, каждая группа будет содержать 
тем больше число процессоров, чем больше шагов интегрирования на базовом шаге 
ей предстоит выполнить для получения аппроксимации решения. 

Из двух рассмотренных вариантов получения экстраполяционной таблицы, имен- 
но второй подходит для пропорционального разбиения процессоров, поскольку не 
требует дополнительной сортировки данных из-за неравномерности их распределе- 
ния. В общем виде проблема разбиения процессорного поля на группы, обеспе- 
чивающие оптимальное решение задачи балансировки загрузки многопроцессорной 
системы, сформулирована в терминах комбинаторной оптимизации. Точное решение 
задач этого класса за приемлемое время невозможно, так как задача является МР- 
сложной, поэтому предлагаются близкие к оптимальному рациональные решения 
приближенными эвристическими методами. 

Сократить время простоя можно при использовании комбинационного способа. 
Идея комбинационного метода состоит в использовании [&/ 7” групп процессоров, 


берется р; = -й; |, а затем оставшиеся процессоры по одному добавляются В 


причем 1-я группа вычисляет 1-ю и (А —1+ /)-ю аппроксимации решения. Для полу- 
чения 1-й аппроксимации необходимо и, раз выполнить обращение к опорному ме- 
тоду решения, для (К—1+1)-й — п;_;./ раз. То есть первая группа должна п; + их 
раз обратиться к вычислению решения по опорному методу, вторая группа — 
п) +пу_| и, наконец, 1 = [%/2] группа — Ик; + Пк 1-1. 

Для разработанных параллельных алгоритмов метода ЛЭР построены схемы 
отображения на параллельные структуры разных топологий: кольцо, тор и гиперкуб. 
Проведенный сравнительный анализ эффективности использования топологий и 
вычислительный эксперимент дают основания утверждать, что менее затратными яв- 
ляются обмены в топологии гиперкуб. 

Определяющим параметром для рассмотренных вычислительных схем методов 
локальной экстраполяции является сложность правой части исходной СОДУ, обес- 
печивающая доминирование вычислений над обменами. Для несложных правых 
частей третья вычислительная схема, объединяющая параллелизм системы и экстра- 
поляции, обладает явным преимуществом, и этот эффект с ростом числа процес- 
соров увеличивается. Для 3ПМО-систем, как и ранее, может быть реализована только 
первая схема, которая имеет низкие показатели ускорения и эффективности. 


Выводы 


Численный эксперимент на базе тестов для СОДУ и проведенный сравнитель- 
ный анализ динамических характеристик параллельных алгоритмов по технологии 
локальной экстраполяции показал, что большая вычислительная сложность после- 
довательных реализаций этих методов и высокая степень внутреннего параллелизма 
являются основанием для эффективного использования их в высокопроизводитель- 
ных мультипроцессорных системах. Апробация параллельных алгоритмов проводи- 
лась на моделях ВС $ПМО, МИМО, кластерной архитектуры с различными топологи- 
ями. В частности, наилучшие характеристики распараллеливания достигаются при 
использовании топологии гиперкуб. Особенно важным является тот факт, что четко 
выявлена область применения параллельных экстраполяционных методов, где они 
практически не имеют конкурентов — высокоточное решение задачи Коши для 
СОДУ со сложными правыми частями. 
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Перспективным направлением дальнейших исследований является применение 
технологии локальной экстраполяции для разработки параллельных методов реше- 
ния жестких задач. Для достижения этой цели необходимо разработать алгоритмы с 
использованием опорных методов, имеющих достаточные характеристики устойчи- 
вости: среди одношаговых методов к таким методам относятся многоточечные блоч- 
ные и полностью неявные методы типа Рунге-Кутты. 
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ТА. Назарова 

Екстраполящин! блоков! однокроков! чисельн! методи розв?язання жорстких задач Коши 
Екстраполящйн! методи розв’язання задач1 Копи для звичайних диференщальних рвнянь мають 
високий стушнь потенщйного паралелзму. Дана стаття присвячена розробщ та анал1зу ефективност! 
паралельних алгоритмв локально! екстраполяци на баз1 блокових явних опорних метод1в. Розроблени 
алгоритми реал1зовано на паралельних системах 13 розподленою пам’яттю 1 тополопею гшеркуб. 
Отримано ошнки часу виконання та обмиев, загальн! накладн! витрати на паралел!зм, прискорення 1 
ефективнисть паралельного розв’язання. 
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